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'nI" I Interstellar dust plays a crucial role in the evolution of the Universe by 

■ assisting the formation of molecules^, by triggering the formation of the 

■ first low-mass stars^, and by absorbing stellar ultraviolet-optical light and 
\ subsequently re-emitting it at infrared/millimetre wavelengths. Dust is 
' thought to be produced predominantly in the envelopes of evolved (age >1 
' Gyr), low- mass stars'^. This picture has, however, recently been brought 
I into question by the discovery of large masses of dust in the host galaxies 

^ ' of quasars^'^ at redshift z > 6, when the age of the Universe was less than 

O ■ 1 Gyr. Theoretical studies^' corroborated by observations of nearby 

r~| ! supernova remnants^ have suggested that supernovae provide a fast 

Qh' and efficient dust formation environment in the early Universe. Here we 

Q ■ report infrared observations of a quasar at redshift 6.2, which are used 

^ , to obtain directly its dust extinction curve. We then show that such a 

^ I curve is in excellent agreement with supernova dust models. This result 

' demonstrates a supernova origin for dust in this high-redshift quasar, from 
which we infer that most of the dust at high redshifts has probably the 
same origin. 



Powerful quasars offer an optimal tool An important class of quasars are the 

for detailed studies of dust in their host Broad Absorption Line (BAL) quasars, 

galaxies out to very high redshifts. Dust whose UV spectrum is characterized by 

extinction is inferred through redden- blueshifted, deep and broad absorption 

ing of the quasar UV/optical continuum features associated with highly ionized 

emission. The extinction curve {Ax = atomic species, which trace powerful out- 

1.086 T\, where tx is the wavelength- flows of dense gas along our line of sight. 

dependent optical depth) of dust associ- Low Ionization BAL (LoBAL) quasars are 

ated with low redshift, mildly obscured characterized by additional low ionization 

quasars has been typically found to be absorption lines, probably associated with 

consistent with that of the Small Magel- higher column densities of gas. LoBAL 

lanic Cloud (SMC)^'^'^^'^^, while for heav- quasars at z < 4 are always significantly 

ily absorbed quasars there are indica- reddened by dust (associated with the out- 

tions that the extinction curve may be flowing gas), and therefore they are ideal 
differentia ' ^6 _ 
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laboratories to investigate the properties 
of dust. 

By means of low resolution near-IR 
spectroscopic observations we have re- 
cently identified a few BAL quasars at 
z 5 - 6 (Maiolino et al. 2004). The 
most distant among them is the LoBAL 
quasar SDSSJ104845.05+463718.3 (here- 
after SDSS1048-F46). At a redshift z = 
6.193 (based on a new medium resolution 
spectrum around MgIIA2798; Maiolino et 
al. in prep), this is the only LoBAL 
at ^ > 5 currently known, and offers a 
unique chance of investigating the dust 
extinction curve at z « 6. This quasar 
was re-observed at higher spectral reso- 
lution with the goal of de-blending the 
deep troughs characterizing its spectrum 
and of estimating the reddening of the 
continuum. Observations were obtained 
in February 2004 with the Near Infrared 
Camera Spectrometer (NICS"*^^) at the 
Telescopio Nazionale Galileo (Canary Is- 
lands). SDSS1048-F46 was observed for 
a total of 3 hours, with a slit of 1.5" 
oriented along the parallactic angle. We 
used a grism covering the spectral region 
0.8 - 1.45/xm. 

Fig. la shows the observed spectrum 
(solid line) smoothed to a lower resolution 
for sake of clarity, and combined with a 
previous low resolution spectrum^''^ of the 
full wavelength range 0.8 — 2.4;um. At 
the quasar's redshift our data cover the 
wavelength range 1200 < Xrest < 3300A 
in the rest frame of the emitted radia- 
tion. The observed spectrum is much 
bluer than LoBAL quasars observed at 
z <4, whose average spectrum is shown 
with the dashed line^^. In particular, 
SDSS1048-I-46 is bluer than any known 
LoBAL quasar^^ at z <4. Since the red 
shape of LoBAL quasars at z<4 is ascribed 
to dust reddening, the much bluer slope of 
SDSS1048 already indicates a large varia- 
tion of dust extinction and reddening from 
z < 4 to 2; Rs 6. 

The continuum at Xrest >1700 A can 



be fitted with a non-BAL (unreddened) 
quasar template with a slope a = —2.1 
{Fx (X A"), which is shown with a dot- 
ted line in Fig.l. Such a blue slope is 
consistent with the intrinsic slope of - 
2.01 inferred by Reichard et al. (2003)^^ 
for BAL quasars. The observed spec- 
trum deviates significantly from the non- 
BAL unreddened template only at Xrest < 
1700 A. While such a deviation was ini- 
tially ascribed to a severe blend of the 
CIV and SilV troughs in the low resolu- 
tion spectrum, the new medium resolu- 
tion spectrum (whose unsmoothed and en- 
larged version is shown in Fig. lb) clearly 
shows that this is not the case. The con- 
tinuum outside the troughs is indeed red- 
der than expected by the extrapolation of 
the longer wavelength spectrum through 
the unreddened quasar template (dotted 
line). If such reddening at short wave- 
lengths is due to dust, then the extinc- 
tion curve must be quite unusual: rela- 
tively fiat at A > 1700 A and steeply rising 
at shorter wavelengths. We can quantita- 
tively derive the extinction curve by using 
the equation Ax = -2.5 log{Fobs/Fintr), 
where Fo6s is the observed spectrum and 
Fintr is the intrinsic spectrum. We avoided 
spectral regions contaminated by emis- 
sion/absorption features by interpolating 
the continuum from the nearby regions. 
The blend of Fell lines^^ makes more un- 
certain the description of the continuum 
in the spectral region between 2300A and 
3050A; however, the latter spectral re- 
gion is not critical within the context of 
this paper, since the most important as- 
pect is that the extinction curve must be 
rather flat in the region between 1700A 
and 3300 A, and increase rapidly at Xrest < 
I700A. For the intrinsic spectrum we used 
the non-BAL template obtained by the 
SLOAN survey-*^^ (this choice is appropri- 
ate since non-BAL quasars at z^-^Q have 
spectra similar to those at z<4)"'^^, with 
slopes ranging from a = —2.1 (the red- 
dest slope compatible with our observed 
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spectrum at Xrest > ITOOA) to a = —2.5 
(only 1% of quasars have slopes bluer than 
this value^^). The resulting extinction 
curve is shown in Fig. 2 (thick solid line 
and shaded region). As expected, the ex- 
tinction curve inferred for this quasar at 
z =6.2 is quite different with respect to 
the SMC curve which applies to quasars 
at z < 4. The dot-dashed line in Fig. lb 
indicates the non-BAL template absorbed 
with the extinction curve inferred by us in 
Fig.2, which nicely matches the observed 
spectrum (except for the emission and ab- 
sorption lines, whose intensity and shape 
depend on the physics of the ionized gas). 
The extinction ^3qqq^ required to match 
the observed spectrum is in the range 0.4- 
0.8 mag. 

The extinction curve inferred from the 
most distant LoBAL can also reproduce 
the shape of the second most distant 
BAL SDSS1044-01, a High Ionization 
BAL (HiBAL)^'^ at z =5.78. The lat- 
ter quasar presents much lower reddening 
than SDSS1048-46 and therefore cannot 
be used to provide tight constraints on the 
extinction curve. Nonetheless, similarly to 
SDSS1048-46, it is characterized a rather 
blue continuum at Xrest > 1700 A and a 
reddening at Xrest < 1700A, which can be 
nicely fitted by using the same extinction 
curve derived above and shown in Fig. 2. 

In order to interpret the observations, 
we have compTitcd the SN dust extinc- 
tion curve by using the model proposed by 
Todini & Ferrara (2001)^ which describes 
dust formation in the ejecta of Type-II 
SNe as a function of the progenitor mass 
and metallicity. In spite of the uncer- 
tainties, the model has been successfully 
applied to interpret the observed prop- 
erties of SN 1987A^ and of the young, 
metal-poor dwarf galaxy SBS 0335-052^^^ . 
A grid of SN models has been consid- 
ered, with different initial stellar progen- 
itor masses (12 < M < 40 M©) 
and metallicities (0 < Z/Zq < if^. 
The resulting SN dust grains are made 



of silicates (Mg2Si04 and MgSiOs), amor- 
phous carbon (AC), magnetite (Fe304) 
and corundum (AI2O3). The typical grain 
sizes range between 5 A to 0.1 /xm, with 
smaller grains being predominantly com- 
posed of silicates and magnetite (sizes ~5- 
15 A) and larger ones by amorphous car- 
bon (sizes «i 100 - 200 A). The pre- 
dicted grain properties for different Type- 
II SNe have then been averaged over the 
stellar Initial Mass Function (IMF). We 
adopted a Salpeter IMF, but even higher- 
mass weighted Larson-like IMFs^^ do not 
change significantly the results ( ^ 2% in 
the extinction curve). 

The SN dust extinction properties have 
been derived using the standard Mie the- 
ory for spherical grains. Grains made 
of AC are the main contributors to ex- 
tinction. We used the optical properties 
for AC produced in an inert atmosphere 
(ACAR'^'^), which is appropriate for the 
SN ejecta environment. Other important 
contributors to extinction were found to 
be Mg2Si04 and Fe304 grains^^'^^. 

Fig.2 shows the resulting extinction 
curves, which are in a much better agree- 
ment with the observations than the usual 
SMC curve. The extinction curves were 
obtained for various initial stellar metal- 
licities. The best agreement is found for 
the Z = 10-2 Zq and 10"^ Zq mod- 
els, although the other models do not dif- 
fer strongly and are still within the ob- 
servational uncertainties. We also show 
the case of dust formed in the ejecta of 
a single 25 Mq, Z = W'^Zq Type-II SN; 
the agreement for this case is excellent at 
most wavelenghts. The plateau between 
1700 A and 3000 A is due to a local min- 
imum between two broad absorption fea- 
tures caused by AC grains, with the added 
flat contribution from Fe304 grains. Sili- 
cate grains contribute to the rise at shorter 
wavelenghts. 

We conclude that our analysis, com- 
bined with observations, provides the first 
direct evidence for dust produced in SN 
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ejecta, rather than in evolved stars, in an 
object at z > 6. In particular, we find 
that dust purely produced by Type-II SNe 
can explain very well the observed extinc- 
tion curve. By implication, we conclude 
that much of the dust seen at high red- 
shifts probably has the same origin. The 
properties of high redshift dust are there- 
fore noticeably different from those found 
at later cosmic times: grains arc typically 
smaller, due to their different formation 
history and to the short time available to 
subsequently accrete heavy atoms and co- 
agulate with other grains. 
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Figure 1. Spectrum of the quasar SDSS1048+46 at redshift 6.2, compared with quasar templates at 
lower redshift. a) The solid line shows the low resolution near-IR spectrum of the Low Ionization Broad 
Absorption Line (LoBAL) quasar SDSS1048+46. It is the composition of a previous low resolution 
A/ A A = 75 spectrum^^ (shown only for A > 1.4/im), and the new medium resolution A/AA = 350 
spectrum smoothed to a resolution of 75 for sake of clarity. Cross-calibration between the two spectra 
is ensured by matching the common parts. Upper labels identify the emission lines while lower labels 
identify the blucshiftod absorption features. The presence of strong absorption from the high ionization 
species CIVA1549 and SIVA1397, along with absorption by the low ionization species A1IIIA1857 and 
MgIIA2798 classifies this as a LoBAL quasar. The missing parts of the spectrum are due to the 
regions of bad atmospheric transmission. The dotted line is the average spectrum of (unreddened) 
non-BAL qxiasars at z < 4 with a slope a — —2.1. The dashed line is the average spectrum of 
LoBAL quasars at 2; < 4. b) Same as above, but where the solid line shows only the (unsmoothed) 
medium resolution spectrum (A/AA = 350) of SDSS1048-I-46. It is important to note the continuum 
level outside the troughs at A < 1700A, which is well below the extrapolation of the continuum at 
longer wavclonghts in the case of no dust reddening (dotted line). Both low and high resolution 
spectra allow the identification of the spectral regions for the continuum fitting (required to derive 
the extinction curve), and more specifically: 1275-1290A, 1325-1340A, 1417-1434A, 1482-1489A, 
I577-I590A, 1567-1677A, 1708-1782A, 1815-1845A, 2015-2270A, 3050-3255A. The dot-dashed line 
in panel b) shows the effect of absorbing the non-BAL template with the extinction curve inferred by 
us and shown in Fig. 2 (thick solid line). 



A SUPERNOVA ORIGIN FOR DUST IN A HIGH-REDSHIFT QUASAR 7 



Kesx [ A ] 

3000 2500 2000 1500 1200 

3.0 



2.5 



2.0 

o 
o 
m 

< 

<^ 1.5 
1.0 



0.5 

3 4 5 6 7 8 

1 / ?i,est [ l^m-' ] 

Figure 2. Extinction curve observed in the quasar SDSS1048+46 a,t z = 6.2 compared with the 
extinction curve observed in quasars at z<4 and with the extinction curve expected from supernova 
(SN) dust. The thick soUd hne shows the extinction curve inferred for SDSS1048+46. The shaded area 
shows the associated uncertainty, which includes the range of intrinsic spectral slopes (this uncertainty 
dominates at short wavelenghts, A < 1700A), the uncertanty on the continuum interpolation (in 
particular around the Fell hump between 2300 A and 3050A), and the noise in the spectrum. The thin 
solid line shows the Small Magellanic Cloud (SMC) extinction curve, which applies to quasars at z < 4 
(including BAL quasars). The other lines indicate theoretical predictions for dust produced by SNe. 
More specifically: dot-dashed, long-dashed, short-dashed and dotted lines represent the extinction 
curves produced by SNe obtained assuming that stars form according to a Salpeter IMF and initial 
metallicities 0, 10~*, 10"'^, and 1 Zq, respectively. Wc also show the case of dust formed in the ejecta 
of a single 25 Mq, Z = 10"* Zq Type-II SN (triple-dot-dashed hne). 




